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INTRODUCTION  WITH  GOALS 

The  goal  of  our  segment  of  the  Unified  Physical  Parameterization  for  Extended  Forecasts  DRI  is  to 
provide  observations  related  to  the  partitioning  of  atmospheric  boundary  layer  (ABL)  fluxes  into  two 
contributions:  eddy  diffusivity  and  vertical  mass  fluxes  (EDMF).  In  addition  to  the  field  campaign, 
the  goal  is  to  validate  the  existing  EDMF  expressions  and  to  propose  alternatives  or  modifications  to 
those  formalisms. 

These  goals  are  being  pursued  using  data  sets  obtained  from  the  CIRPAS  Twin  Otter  flown  in  two 
differing  regions:  off  the  coast  of  California  near  Monterey  and  over  portions  of  the  Dugway  Proving 
grounds  in  Utah. 

OBJECTIVES 

The  general  objective  is  to  collaborate  with  NRL  (Wang)  and  UCI  (Khelif)  to  use  airborne  sensors  to 
probe  the  ABL  for  organizing  structures  such  as  OLEs  and  LLJs  and,  using  models  when  appropriate, 
construct  algorithms  that  would  relate  space-based  sensor  observables  to  the  varying  structure 
dependent  ABL  fluxes.  Two  regions  have  been  used  to  collect  data  for  the  evaluation  and  modification 
of  the  EDMF  formalisms: 

1.  In  areas  of  neutral  to  stable  ABLs  capped  by  cloud  cover  ranging  from  zero  to  100%,  measure 
the  fluxes  at  the  air-sea  interface  using  the  TODWL  and  CTV  instruments  flown  together  near 
the  California  coastline  in  September  2012. 

2.  Using  the  MATERHORN  experiment  at  Dugway,  UT  in  October  2012  as  an  opportunity  to 
expand  the  DWL  data  base  to  include  ABL  surveys  in  complex  terrain  as  a  comparison  to  the 
ABL  over  cold  water.  Note  that  the  MATERHORN  was  co-funded  by  ONR  (Ferek). 

As  this  research  has  progressed,  we  have  added  a  modeling  component  to  better  understand  the 
implications  of  ABL  structures  being  revealed  by  the  TODWL  during  the  September  and  October  2012 
field  campaigns. 

APPROACH 

The  work  on  objective  one  above  continues  as  we  examine  several  cases  in  2013  where  the  CTV  was 
operating  adequately  and  the  ABL  contained  OLEs  and  LLJs  needed  for  our  investigation.  The  first 
important  step  was  to  navigate  the  data  from  the  two  systems  so  that  their  data  sets  could  be  joined. 
This  was  not  as  easy  as  it  seemed  since  the  CTV  did  not  fly  faithfully  directly  below  the  Twin  Otter. 

The  second  step  is  to  compute  heat,  moisture  and  momentum  fluxes  for  both  the  CTV  and  Twin  Otter 
flux  sensors.  The  computation  of  fluxes  from  the  CTV  will  be  done  by  both  UCI  (Khelif)  and  SWA 
(Emmitt). 

The  third  step  is  to  determine  the  correlations  between  ABL  structures  and  the  fluxes.  The  end  product 
is  envisioned  as  being  a  set  of  PDFs  for  mass  flux  contributions  as  a  function  of  regional  model 
representations  of  the  BL  including  space-based  observations  that  can  be  used  to  infer  organized 
structures  such  as  LLJs  (ASCAT  and  CALIPSO)  and  OLEs  (MODIS,  VIIRS  and  GEOS). 
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Foster  has  begun  modifying  his  OLE  model  to  better  understand  the  implication  of  “stacked  rolls” 
being  detected  by  the  TODWL  during  flights  just  off  the  California  coast.  Emmitt  and  de  Wekker  are 
using  the  WRF  and  COAMPs  models  to  test  out  their  sensitivities  to  changes  in  the  EDMF  related  to 
our  field  data.  We  will  use  the  SCM  provided  by  NRL  if  appropriate. 

To  meet  the  second  objective,  we  processed  several  thousand  wind  profiles  obtained  with  TODWL 
over  Granite  Mountain  at  Dugway  Proving  grounds,  UT.  The  primary  data  products  derived  are  the 
three  wind  components,  turbulence  on  the  scales  of  the  illuminated  sample  volumes  and  aerosol 
structures  (PBL  depths,  OLEs,  incipient  clouds,  etc.) 

After  processing  the  TODWL  data  for  7  MATERHORN  missions,  comparisons  of  u,v  and  vertical 
velocities  are  compared  to  those  predicted  by  a  WRF  model. 

A  final  effort  will  be  made  to  relate  the  fluxes  measured  at  flux  tower  sites  with  those  measured 
onboard  the  Twin  Otter  during  flybys.  During  those  passes  the  winds  and  aerosols  were  mapped  by  the 
TODWL.  These  data  sets  will  be  explored  for  relationships  between  the  BL  organization  and  flux 
intensities. 

WORK  COMPLETED 

In  the  final  half  of  >ear  2  on  this  project,  we  have  processed  more  than  50  hours  of  TODWL  data  from 
both  the  Monterey  and  DPG  areas.  We  have  begun  combining  the  information  from  the  TODWL,  CTV 
and  Twin  Otter  sensors  to  establish  the  relationship  between  the  local  fluxes  and  the  energetic  of  LLJs 
and  OLEs.  We  have  also  begun  to  modify  the  EDMF  (Eddy  Diffusivity  and  Mass  Flux) 
parameterization  to  account  for  the  differences  between  thermally  driven  convection  and  dynamically 
driven  vertical  transports. 

The  WRF  model  has  been  setup  and  run  for  the  7  MATERHORN  missions  and  comparisons  between 
the  TODWL  wind  profiles  and  the  model  profiles  have  been  completed  for  one  day. 

The  OLE  model  has  been  modified  to  investigate  “stacked  OLEs”  seen  in  the  April  2007  datasets. 

RESULTS 

Most  of  the  results  of  our  last  year’s  research  effort  were  reported  at  the  DRI  Workshop  in  Monterey  in 
August  2013  and  at  a  Coherent  Laser  Radar  Conference  in  June  2013.  Both  presentations  are 
appended  to  this  report. 

The  two  presentations  are  entitled: 

PBL  structures  governing  the  Unified  Physical  Parameterization 
formulation  (Emmitt,  Foster  and  deWekker) 

Airborne  DWL  investigations  of  flow  over  complex  terrain 
(MATERHORN  2012)  (Emmitt,  Godwin,  Greco  and  de  Wekker) 
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In  addition  to  the  presentations  provided  as  appendices,  a  report  compiled  by  de  Wekker  is  also 
attached.  The  material  presented  by  de  Wekker  combines  MATERHORN  research  results  that  were 
funded  by  ONR.  ARO  and  NSF.  This  work  lays  the  foundation  for  UPP  related  investigations  that  will 
be  carried  out  in  Year  3. 

RELATED  PROJECTS 

ONR  contract  to  study  the  utilization  of  Doppler  wind  lidar  (DWL)  data  to  quantify  the  contribution  of 
organized  large  eddies  (OLEs)  to  fluxes  in  the  marine  boundary  layer  (DYNAMO  DRI  (S.  Harper, 
TPOC). 

PLANS  FOR  YEAR  3 

We  have  four  primary  goals  for  year  3  of  this  effort. 

1 .  Finish  the  development  of  a  methodology  for  joining  airborne  DWL  data  and  in  situ  flux  data 
to  achieve  an  understanding  of  how  the  organized  ABL  structures  modulate  ABL  fluxes. 

2.  Apply  that  methodology  to  the  DWL  data  sets  already  identified  and  in  some  cases  analyzed. 

3.  Develop  an  expression  for  mass  fluxes  by  OLEs  in  the  MBL  that  is  suitable  for  use  in 
numerical  model  parameterizations.  This  will  be  derived  from  flux  data  measured  directly  by 
the  CTV  (Monterey)  and  towers  (Dugway)  and  associated  with  organized  ABL,  structures 
mapped  from  the  DWL. 

4.  Use  the  COAMPS,  WRF  and  SCM  to  investigate  the  sensitivity  of  these  models  to  the  ranges 
of  v  alues  for  the  “MF"  term  as  inferred  from  the  Monterey  and  Dugway  data  sets. 

We  will  pursue  publishing  our  results  should  the  findings  of  this  research  merit  this  effort. 

We  will  also  prepare  a  work  plan  and  budget  for  years  4  and  5.  This  plan  will  be  submitted  by  the  end 
of  2013.  Fhe  primary  issues  to  be  considered  will  be: 

Can  we  justify  a  second  field  project  in  FY  2014  that  will  focus  upon  the  lidar  resolved  MBL 
structures  and  their  impacts  on  fluxes?  In  2012  the  lidar  was  flown  on  just  one  good  mission 
due  to  weight  considerations  and  the  desire  for  5hr  flights.  This  would  not  be  the  case  in  a 
proposed  second  series  of  TODWL/CTV  flights. 

What  space-based  data  sets  can  be  used  to  classify  the  type  and  degree  of  “organized”  BL 
structures  and  thereby  adjust  the  coefficients  in  the  EDMF  being  used  in  global  models? 

What  is  the  best  flux  parameterization  that  recognizes  the  differences  in  the  mass  flux  traced  to 
dr}'  thermal  convection  vs.  dynamically  driven  BI.  circulations? 
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Airborne  DWL  Investigations  of  Flow  over  Complex  Terrain 
(MATERHORN  2012) 

Paper  presented  at  the  Coherent  Laser  Radar  Conference,  June  2013 
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flow  over  complex  terrain 
(MATERHORN  2012) 
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(University  of  Virginia) 


CLRC  2013 
Barcelona ,  Spain 
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Outline 


Science  objectives  of  the  MATERHORN 
Experiment  summary  (TODWL  data) 
WRF  model 

DWL  soundings  on  Google  Earth 
WRF  and  TODWL  data  comparisons 

Funding,  provided  by: 

-  Office  ofNavaj  Research  (Fereki 

-  Army  Research  Office  (Videen) 
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*  twin  Otter  Doppler  Wind  Lidar 
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TWOLF 


Tornado-chasing  Wind  Observing  Lidar  Facility 


Pallet  Wind  Observing  Lidar  Facility  (PWOLF) 


Attribute 

Performance  Metric 

Comments 

LOS  resolution  (applies 
to  vertical  profiles  of 

3D  winds  as  well) 

50  m 

Range  resolution  to 
hard  targets  (ground  or 
dense  cloud)  can  be 
better  than  10  meters. 

U,V,W  resolution 

<  1 0  cm/s 

<  5  cm/s  for  stationary 
groundbased  operations 

Maximum  range 

6  -30  km 

Very  dependent  upon 
aerosols 

Time  to  complete  full 
conical  scan  for  wind 
profiles 

~  20  sec 

Sampling  frequency 

100  Hz 

Integration  of  several 
shots  is  typical  to 
improve  range 
performance 

CIRPAS  Twin  Otter  at  MATERHORN 

MATERHORN  is  a  multi-agency  research  project 
to  better  understand  and  model  atmospheric 
circulations  in  complex  terrain. 

—  Provide  3D  context  for  complex  terrain  flow  studies 
being  conducted  with  towers,  rawinsondes  and  UAVs, 

The  MATERHORN  also  provided  the  opportunity 
to  pursue  Unified  Physical  Parameterization 
objectives  with  evaluation  of  EDMF 
parameterization  over  land  in  the  presence  of 
organized  structures  in  daylight  and  nocturnal 
flows  and  katabatic  flows. 


Flight  profile 


Navy  Twin  Otter  based  out  of  Salt  Lake  City 

—  ~  20  minute  to  Granite  Mountain 

—  climb  to  12K  feet  (~5K  feet  above  peaks) 

Twin  Otter  in  Utah  between  5  October  and  18 
October,  2012 

Missions  lasted^  4  hours 

—  7  missions  yielded  ~3200  wind  profiles  between 
surface  and  3400  meters. 


TOD.WLd.ata  products 


Downward  conical  scans  (12  point  step  stare) 

Off-nadir  angle  of  20  degrees 

—  20  -25  seconds  for  full  360  scan  (  ~  1  -1.2km) 

—  U,V,VV  with  50  m  vertical  resolution 

•  precision  (U,V)  ~  .10  m/s ;  precision(W)  ~  .15  m/s 

—  SNR  (aerosols) 

Nadir  samples 

—  5  seconds  between  conical  scans 

—  50  m  vertical  resolution  with  w  precision  ~  .10 

Structure  prospecting 

—  Straight  ahead  and  down  3  -6  degrees 


3-D  Doppler  Wind  Profiles  near 
Dugway,  UT  during 
MATERHORN  2012 
10/06/2012 
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MATERHORN  2012  TODWL  N-S  Legs  on  October  9,  2012 

1601  - 1641 


*  -  Provided  by  Simpson  Weather  Associates,  NRL,  ONR,  ARL  and  ARO 


MATERHORN  2012  TODWL  E-W  Legs  on  October  9,  2012 
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*  -  Provided  by  Simpson  Weather  Associates,  NRL,  ONR,  ARL  and  ARO 


Vertical  wind  component 
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Mean  VLOS  by  Heading/WD 
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VLOS  by  Horizontal  Wind  Speed 
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Wnadjr  as  function  of  roll  and  pitch 
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WRF  Model 

•  WRF*  Version  3.3  to  generate  3D  winds 

•  3  Nested  Grids  -  9  km,  3  km,  1  km 

•  Dugway,  Utah  area 

—  Used  for  previous  DoD  studies, 
flight  test  campaigns 

—  Reasonably  representative  of  theaters 
of  interest 

•  Domain:  35  km  x  35  km  horizontal  x  42  vertical  levels 

-  10  minute  time  steps  available  and 
utilized  for  some  of  the  analysis 


*  Weather  Research  Forecasting  Model 


N/S  legs  flown  at  3400m  (MSL)  on  10/9/12 
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Wind  Speed  West  of  Granite  Mountain  -  10/09/2240Z 
Taken  from  WRF  MODEL  1  km  Domain 
Matching  N-S  Lidar  Leg  of  2230-2240 
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Wind  Direction  West  of  Granite  Mountain  -  10/09/2238-40Z 
TODWL  Lidar 
Elevation  -  1300  m  amsl 


Height  (m  AGL) 


Wind  Speed  on  East  Slope  of  Granite  Mountain  -  10/09/2210Z 
Taken  from  WRF  MODEL  1  km  Domain 
Matching  N-S  Lidar  Leg  of  2213-2216 
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Wind  Speed  on  East  Slope  of  Granite  Mountain  - 10/09/221 3-1 5Z 

TODWL  Lidar 
Elevation  ~  1315m  amsl 
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Wind  Direction  on  East  Slope  of  Granite  Mountain  - 10/09/221 OZ 

Taken  from  WRF  MODEL  1  km  Domain 
Matching  N-S  Lidar  Leg  of  2210 
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Wind  Direction  on  East  Slope  of  Granite  Mountain  - 10/09/221 3-1 5Z 

TODWL  Lidar 
Elevation  ~  1315m  amsl 


Height  (m  AGL) 


Wind  Speed  East  of  Granite  Mountain  - 1 0/09/221 OZ 
Taken  from  WRF  MODEL  1  km  Domain 
Matching  N-S  Lidar  Leg  of  2213-2216 
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Wind  Speed  East  of  Granite  Mountain  -  10/09/2206-08Z 
TODWL  Lidar 
Elevation  ~  1316m  amsl 
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Wind  Direction  East  of  Granite  Mountain  - 10/09/221 OZ 
Taken  from  WRF  MODEL  1  km  Domain 
Matching  N-S  Lidar  Leg  of  2213-2216 
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Summary 

Weather  patterns  provided  several  different  flow  patterns 
over  study  area  (Dugway  Proving  Grounds) 

TODWL  performed  very  well  with  no  unscheduled  down 
times. 

Other  data  were  provided  by  UAVs,  tethered  balloon,  ground- 
based  DWLs,  instrumented  towers,  rawinsondes  and  smoke 
releases. 

TODWL  data  is  supporting  model  (WRF/NCAR  and  WRF/Pu) 
validation  and  development  being  conducted  at  UVa  and 
SWA. 

TODWL  data  combined  with  flux  tower  data  addressing  KDMF 
parameterization  schemes 

Data  also  being  used  to  improve  algorithms  for  Precision 
AirDrops. 
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PBL  Structures  Governing  the  Unified  Physical  Parameterization  (UPP)  Formulation 
Paper  presented  at  the  project  workshop  in  Monterey,  CA,  August  2013 


PBL  structures  governing  the 
Unified  Physical  Parameterization 

formulation 

G.  D  Emmitt,  SWA 
R.  Foster,  API  LJWash 
S,  De  Wekker,  UVA 


Overview 


Objectives  of  the  September  2012  field  campaign 
Description  of  observing  system  and  data 
Summary  of  DWL  data  collection  out  of  Monterey,  CA 
Summary  of  DWL  data  collection  at  DPG,  Utah  in  October 
Best  day  data  09/30/12  as  example 

—  Flight  path 

—  Clouds 

—  Vertical  position  record 

•  Twin  Otter  and  CTV 

—  Vertical  soundings 
Additional  case  studies 

Modeling  of  rolls  and  their  relative  contributions  to  EDMF 
Near  term  plans 


Objectives 


Extend  prior  investigations  (2001-2008)of  LLJs 
and  OLEs  in  the  MBL  and  PBLs. 

Investigate  and  characterize  the  presence  of 
rolls  (OLEs)  at  the  boundaries  of 
stratocumulus  topped  MBLs. 

Study  the  potential  impact  on  the  develoment 
and  implementation  of  the  EDMF  into  forecast 
models. 


The  observing  systems  and  strategies 


TODWL 

Twin  Otter  Doppler  Wind  Lidar 


A  SWA 


Attribute 

Performance  Metric 

Comments 

LOS  resolution  (applies 
to  vertical  profiles  of 

3D  winds  as  well) 

50  m 

Range  resolution  to 
hard  targets  (ground  or 
dense  cloud)  can  be 
better  than  10  meters. 

U,V,W  resolution 

<  1 0  cm/s 

<  5  cm/s  for  stationary 
groundbased  operations 

Maximum  range 

6  -30  km 

Very  dependent  upon 
aerosols 

Time  to  complete  full 
step  stare  conical  scan 
for  wind  profiles 

~  20  sec 

12  point  step  stare  with 
.5  -2  second  dwells 

Sampling  frequency 

100  Hz 

Integration  of  several 
shots  is  typical  to 
improve  range 
performance 

TODWL  under  flight  of  WINDSAT  on  18  April  2007  6:07  pm  PST 

Monterey  Bay,  California 
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Some  additional  information 


Feature  prospecting  uses  a  very  shallow  angle 
below  the  horizon  (~  -1  -3  degrees  for  a  300m 
flight  altitude). 

—  Results  in  ~  2  m  vertical  resolution  and  50  m 
horizontal  resolution  with  ~10  meter  sliding 
sample. 

—  It  takes  ~  40  seconds  to  profile  100  meters  below 
the  aircraft. 


Processing  lidar  returns  in  the  spectral 

domain 


Data  summary  for  September/October 
UPP  field  campaign  with  TODWL 

•  Three  quality  flights  with  the  DWL 

•  One  quality  co-flight  with  the  CTV  (9/30). 

•  Two  5  hour  ferry  flights  between  Monterey,  CA 
and  Salt  Lake  City,  Utah. 

•  Seven  4.5  hour  flights  over  DPG,  Utah  during  the 
MATERHORN 

-  Co-funded  by  ONR  and  ARO 

-  3500  vertical  profiles  used  to  validate  WRF  models 

-  14  low  level  flights  dedicated  to  prospecting  for  OLEs 
in  the  vicinity  of  flux  towers. 


Pallet  Wind  Observing  Udar  Facility  (PWOLF) 


(s/wj)  poods  pu|m 


Summary  of  the  7  research  flights  conducted  during  MATERHORN  in 
support  of  the  UPP  and  MATERHORN  projects. 

Winds  at  700  and  850  mb  provide  an  indication  of  the  magnitude  of  the 
synoptic  forcing  during  the  various  flight  days. 

Flights  were  conducted  during  MATERHORN  IOP's  4,  5,  6,  and  7. 


A  major  effort  is  underway  to  investigate  the  presence  of 
organized  convective  structures  under  these  different 
conditions  and  the  contribution  of  these  structures  to  vertical 
heat  and  momentum  fluxes.  These  investigations  include  a 
detailed  analysis  of  surface  fluxes  obtained  by  in-situ 
turbulence  sensors  on  the  many  towers  that  were  deployed 
during  MATERHORN.  A  post-doc  (Sandip  Pal)  and  a  graduate 
student  (Mark  Sghiatti)  at  UVA  will  contribute  to  this  effort 
during  the  next  year. 


9/30/12  UPP:  Monterey  case  study 


9/30/12  case  study 

•  Process  lidar  data  in  search  of  organized 
aerosol/wind  structures  below  the  Twin  Otter 
flight  path 

•  Process  Twin  Otter  "CABIN"  data  for  time  series 
of  u,  v,  w,  q,  and  0. 

•  Process  CTV  data  for  u,v,w,q  and  9. 

•  Match  up  times  and  then  features  from  the 
TODWL  and  CABIN  data  sets  near  flight  level. 

•  Match  up  times  and  features  from  TODWL  and 
CTV  at  CTV  cruise  levels. 


Flight  path 


Time  (LST)  and  Location  of  CTV/Lidar  Data  and  1045  LST  GOESW  Imagery 


u  jis  (L  r)  and  Location  ci  CTV/Lidar  Data  and  1145  LST  GOESW  Imagery 
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Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 
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Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Twin  Otter:  30  Sep-2012  16:20:39  to  19:15:59  UTC 


Twin  Otter:  Racetrack 


Twin  Otter:  Racetrack 


LOS  Velocity  093012  100^09  Az.  Range  -8  to  2  El.  Range: -3  to  0 
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LOS  Velocity  093012  110626  Az.  Range.  9  to  1 1  El.  Range. -3  to -1 


Spectral  Width 
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LOS  Velocity  09301 2  1 1 1 1 26  Az.  Range.  8  to  1 1  El.  Range.  -3  to  -2 
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LOS  Velocity  093012  112126  Az.  Range. -6  to -4  El.  Range:  -4  to  -2 


093012-102845  Heading:  98.1174 


Longitude 


Height  (m  MSI) 


LOS  Velocity  093012  102845  Az.  Range.  -7  to  -6  El.  Range: -3  to -2 


Spectral  Width 
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Summary  of  segment  statistics 


Flight 
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Heading 
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Is  Double-OLE  Interpretation 

Reasonable? 

Simplify  Problem: 

—  Neutral  stratification 

•  Shear  effects  only 

•  Omit  effect  of  stratification  above  jet 

•  No  surface  buoyancy  flux 

—  Elevated,  thin  baroclinic  layer 

—  Variable  K(z)  (akin  to  MRF/YSU) 

Non-linear  stability  model 
Interacting  triads 


Twin  Otter:  30-Sep-2012  16:20:39  to  19:15:59  UTC 


Mean  wind  at  flight  level  is  generally  from  the  north 

Twin  Otter:  Racetrack  Twin  Otter.  Racetrack 


TO  profiles  in  "racetrack" 


15 


Z  (ND) 


Simplified  Profile  (thin  baroclinic  layer  above) 


Wavelength:  Larger-scale  structures  ~  1.5  to  2.5  km 
Smaller-scale  structures  ~  300  to  700  m 

Along-roll  velocity  Perturbations:  +/-  20-25  % 

Note;  largest  near  surface 


Orientation:  Typically  aligned  along  shear 
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Dominant  mode 


90  2 


120 


240 


300 


270 


Secondary  mode  associated  with 
above-jet  layer 


90  2 


Single-Wave  Roll  Theory 
Nonlinear  Stability 


"Stretch"  eigenvalue,  A0,  in  powers  of  nonlinear  amplitude,  Aft). 
Expand  eigenfunction.  q1Q,  in  harmonics  of  fundamental  wavenumber, 
a,  and  forced  modifications 

•  Forced  fundamental  modifications  are  orthogonal  to  linear  mode 

•  Determine  Landau  Coefficients  (the  A.;) 

Estimate  equilibrium  Amplitude  (dA/dt  =  0)  &  structure,  q  =  [u,v,w,Tr 


\\j,  a  =  1 .6, 


Primary  Mode 


8  =  20° 


•  Mostly  in  the  sub-jet  layer 

•  Along-roll  velocity 
strongest  near  jet  as 
compared  to  std.  rolls 

•  Oriented  at  large  angle  to 
dominant  rolls 


WW,  a  =  1 .6,  s  =  20° 


a  4 


0  5~  10  15  20  25 

X  (ND) 


Weaker,  secondary  mode 


•  Similar  to  STD  rolls,  but 
based  on  jet  top 

•  Along-roll  velocity  max 
near  jet 

•  Oriented  at  large  angle  to 
dominant  rolls 


\\ j,  a  =  0.5,  s  =  80° 
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0 - 1 - 1 - 1 - 1 - - 1 

0  5  10  15  20  25 

X  (ND) 


Resonant  Triad  Interaction  Between  OLE 

•  a  =  (3  +  y  (mode  A,  mode  B,  mode  C) 
wavenumber  relationship 

•  Require  at  least  one  wavenumber  at  fastest 
growing  mode  (A) 

•  Exchange  of  energy  between  OLE  modes 

•  Often  energizes  weak  modes 

•  At  present,  require  co-linear  modes  (can  be 
relaxed) 
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Truncated  Model 

Amplitude  (real)  and  Phase  (imaginary) 


1  dA 
A  dt 
1  dB 


B  dt 
1  d£ 
C  dt 


.  ddA 

i — a 
dt 

,  ddD 
l 


=  ao 


dt 
.  <20 


bn 


l 


dt 


c  — 


=  c0 


+  [a2A2  +  a2B-  +  a4C2] 
+  [b ,  A2  +  b3B2  +  b4C2] 
+  [c?A  +  c3B  +  c4C  ! 


4>  =  0^  —  —  0r  (Wave  phase  imbalance) 


a  =  (3  +  y  (resonant  triad  wavenumbers) 


The  aj,  bj,  q  are  complex,  generalized  Landau  coefficients,  calculated  via  an 
orthogonalization  assumption  (nonlinear  wave-wave  &  wave-mean  flow 
interactions) 

Highest-order  (bracketed)  terms  force  equilibrium;  dominated  by  single-wave 
contributions  (a2A2,  b3B2,  c4C2) 

Lower-order  phase  coupling  allows  inter-scale  energy  transfer,  enhances  growth 

RATE  OF  SLOWEST-GROWING  MODE,  ESPECIALLY  DURING  QUASI-LINEAR  PHASE 


ND  Amplitude 


How  is  the  slowly-growing  mode  energized? 
Examine  quasi-linear,  early-time  behavior 


Accelerated  growth  of  slow  mode  (C) 


Starts  when  (c|>  -  4>cl )  2ttn 


Phase  Imbalance: 
4>  =  0A  +  ®c 


QL  Landau  Coefficient: 


i<t>y 
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Z  (ND)  Z  (ND) 


Xa  =  5.03  (ND);  X  =  7.85  (ND);  X  =13.96  (ND) 
W  (shading)  and  vy  (contour),  time  =  0150  (ND) 


U1  (shading)  and  \\i  (contour) 


x  (ND) 


Sub-Summary 

•  It  appears  that  a  PBL  jet  can  induce  "stacked"  OLE 

—  Caveat:  test  case  is  neutral  stratification 
•  In  model,  but  examine  neutral  first 

•  CO-LINEAR  OLE  TRIADS  CAN  GENERATE  COMPLEX 
STRUCTURE 

•  There  is  a  weaker,  larger-scale,  oblique  above-jet 

MODE 

-  Can  it  be  energized  by  nonlinear  triad  mechanism? 

•  Future  work 

-  Modify  triad  code  for  non-co-linear  OLE 

—  Include  stratification  and  thermal  wind  effects  on  OLE 


Revise  and  Add  to  Overall  Summary 


•  EDMF-type  schemes  attempt  to  capture  non-local 
contributions  to  PBL  fluxes 

-  Implicitly  assumes  such  transport  is  due  to  narrow, 
skewed  updrafts 

—  Quite  unstable  stratification 

•  OLE  rolls  are  very  common 

-  low  skewness 

-  Likely  significant  flux  contributions 

—  Near-neutral  to  moderately  unstable  stratificaton 

•  Possible  also  slightly  stable? 
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Extra  Slides 


Truncated  3-Mode  Roll  Solutions 


qa  = 


% 


qv  = 


:  contributions  from  single-wave  theory;  e.g.  q2  a  =  qn  a . 

•  BLUE:  new  wave-wave  &  wave-mean  flow  interaction  contributions. 

•  Red:  Low-order  phase-coupling  terms. 

•  Also:  mean-flow  modifications  due  to  each  wave. 
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Standard  Non-Linear  Single-Wave  PBL  Roll  Model 


Table  5,1  Contributions  to  the  nonlinear  perturbation  up  to  the  fifth  Landau 
Coefficient 


Truncated  Contributions 
to  Multi-Wave  Roll  Model 


q=  [u,  v,  w,  T] 


T 


0  +  A2q01 

Aq10  +  0  +  A3qlA 

0  +  A2q20 

0  +0 


(mean  flow  modification) 
(fundamental  wavelength) 
(1st  harmonic) 

(2nd  harmonic) 


+  A3q30 
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Low-Order  Truncation  Errors 


Re  =  500.  alpha  =  0.500.  epsilon  =  16.00  deg 


Re  =  1 50.  alpha  =  0.540.  epsilon  =  1 1 .00  deg 


Mean-flow  Modification 


Amplitude  Estimation 


Low-order  truncation  problems: 

•  Over-estimated  amplitude 

•  "S-shaped"  MF  modification 
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Future  plans 
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Attachment  3 

Report  of  University  of  Virginia  Project  Efforts 


Investigation  of  the  representation  of  OLEs  and  terrain  effects  within  the  coastal  zone  in 
the  EDMF  parameterization  scheme:  an  airborne  Doppler  wind  lidar  perspective. 

Stephan  F.J.  De  Wekker 
University  of  Virginia 

Reporting  period:  1  July  2012-30  June  2013 

The  planned  tasks  for  year  2  of  this  pro  ject  that  involved  subcontractor  UVA  were: 

Task  2.1  -  Planning  of  flights  for  field  program  (SWA.  UW.  and  UVA)  to  optimize 
evaluation  of  EDMF  parameterization. 

fask  2.2  -  Field  program  participation.  UVA  will  provide  radiosonde  launches. 

Task  3.1  -  Processing  and  analysis  (including  algorithm  and  software  development)  and 
of  new  TODWL  data  after  completion  of  field  program 

Task  3.7  -  Analyze  TODWL  data  and  compare  with  Radiosonde  date  (UVA  and  SWA) 
Accomplishments: 

During  year  2.  a  majority  of  the  effort  was  in  planning  and  conducting  of  TODWL  research 
flights.  The  UVA  contribution  was  focused  on  investigating  OLEs  and  their  representation  in  the 
EDMF  scheme  over  land.  Flights  were  conducted  for  this  part  over  the  Salinas  Valley  (1  research 
flight)  and  7  flights  were  conducted  over  Dugway  proving  Ground  in  Utah  as  part  of  the  ONR 
funded  MATERHORN  field  campaign.  UVA  provided  and  launched  about  20  radiosondes  for 
these  field  campaign.  The  TODWL  flight  pattern  that  was  developed  for  Dugway  Proving 
Ground  consisted  of  straight  north-south  and  east-west  legs  at  12.000  ft  and  a  sequence  of  low 
level  flights.  The  legs  were  centered  on  Granite  Mountain,  a  mountain  ridge  of  about  1000  m 
high.  A  typical  flight  pattern  is  shown  in  Fig.  1 .  During  the  north-south  and  east-west  legs,  the 
airborne  lidar  was  scanning  in  a  step-stare  pattern  obtaining  12  radial  velocity  profiles  during  one 
full  rotation.  These  data  were  processed  to  obtain  u.  v.  and  w  wind  components  at  50-m  vertical 
resolution.  Between  tw  o  full  rotations  of  the  scanner,  the  lidar  was  pointing  straight  down  (nadir 
stares)  resulting  in  accurate  measurements  of  the  vertical  velocity.  An  example  of  the  vertical 
wind  profiles  during  a  north-south  leg  is  shown  in  Fig.  2.  During  this  particular  flight,  upper- 
level  winds  were  generally  from  southwesterly  directions  between  5  and  10  m/s.  Around  and 
below  the  height  of  Granite  Mountain  at  2  km.  a  drastic  change  in  wind  speed  and  direction  can 
be  observed. 
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1500  2000  2500  3000  3500  4000 


Longitude  (degr.) 


Figure  1:  Flight  pattern  during  RF04  with  the  colors  indicating  the  height  of  the  aircraft.  Typically,  a  north-south  and  east  west 
flight  pattern  was  conducted  to  obtain  vertical  profiles  of  wind  speed  and  direction,  followed  by  low  level  flights  for  the 
prospecting  of  organized  convective  structures. 
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Figure  2:  vertical  profiles  of  horizontal  w  ind  along  the  easternmost  N-S  leg  shown  in  Fig.  1. 


The  low  level  flights  were  designed  to  prospect  for  organized  convective  structures  and  for 
quantifying  vertical  heat  and  momentum  fluxes  associated  with  these  structures. 


Steve  Greco  (SWA)  and  Stephan  De  Wekker  (UVA)  have  been  performing  simulations  on  a  64 
processor  cluster  at  UVA.  In  addition.  De  Wekker  has  started  a  collaboration  with  the  Research 
Application  Laboratory  at  the  National  Center  for  Atmospheric  Research  (NCAR-RAL)  to 
perform  ‘very  large  eddy  simulations'  (VLES)  with  a  300  m  grid  spacing  for  the  various  days  of 
the  research  flights.  An  initial  comparison  of  the  winds  in  a  horizontal  cross  section  shows  that 
some  important  features  of  the  observed  winds  are  captured  by  the  VLES  (Fig.  3).  A  more 
detailed  analysis  and  evaluation  the  VLES  and  of  a  mesocale  model  in  which  an  EDMF  scheme 
is  implemented  is  planned  for  the  third  year  of  this  project.  Additionally,  we  are  planning 
idealized  numerical  simulations  using  WRF-LES  to  simulate  OLE?s  (in  particular  rolls)  that  will 
be  used  to  test  a  theoretical  model  for  OLE  dynamics  that  has  been  developed  by  collaborator 
Foster  (UW). 
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Figure  3:  comparison  between  observed  (left)  and  simulated  (right)  winds  at  1800  m  AGL  during  RF03. 


A  summary  of  the  flights  that  were  conducted  in  support  of  UPP  and  MATERHORN  is  shown  in 
Figure  4.  Six  out  of  the  seven  flights  were  conducted  under  quiescent  to  moderate  condition  as 
indicated  by  the  700  and  850  mb  winds  on  the  flight  days.  RF07  was  conducted  under  stronger 
synoptic  forcing  than  other  days.  A  major  effort  w  ill  be  put  in  the  third  year  of  this  project  to 
investigate  the  presence  of  organized  convective  structures  under  these  different  conditions  and 
the  contribution  of  these  structures  to  vertical  heat  and  momentum  fluxes.  These  investigations 
w  ill  also  include  a  detailed  analysis  of  surface  fluxes  obtained  by  in-situ  turbulence  sensors  on 
the  many  towers  that  were  deployed  during  MATERHORN.  A  post-doc  (Sandip  Pal)  and  a 
graduate  student  (Mark  Sghiatti)  at  UVA  will  contribute  to  this  effort  during  the  next  year. 


Figure  4:  Summary  of  the  7  research  flights  conducted  during  MATERHORN  in  support  of  the  UPP  and  MATERHORN 
projects.  Winds  at  700  and  850  mb  provide  an  indication  of  the  magnitude  of  the  synoptic  forcing  during  the  various  flight  days. 
Flights  were  conducted  during  MATEREIORN  lOP's  4.  5.  6,  and  7. 


Publications  related  to  the  project: 
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